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Abstract. Effect of the ring current ions in the real part of electromagnetic ion 
cyclotron wave dispersion relation is studied on global scale. Recent Cluster observations 
by Engebretson et al. [2007] showed that although the temperature anisotropy of 
energetic (> 10 keV) ring current protons was high during the entire 22 November 2003 
perigee pass, electromagnetic ion cyclotron waves were observed only in conjunction with 
intensification of the ion fluxes below 1 keV by over an order of magnitude. To study the 
effect of the ring current ions on the wave dispersive properties and the corresponding 
global wave redistribution, we use a self-consistent model of interacting ring current 
and electromagnetic ion cyclotron waves [ Khazanov et al., 2006], and simulate the 
May 1998 storm. The main findings of our simulation can be summarized as follows: 
First, the plasma density enhancement in the night MLT sector during the main and 
recovery storm phases is mostly caused by injection of suprathermal plasma sheet H + 
(^: 1 keV), which dominate the thermal plasma density. Second, during the recovery 
storm phases, the ring current modification of the wave dispersion relation leads to a 
qualitative change of the wave patterns in the postmidnight-dawn sector for L > 4.75. 
This “new” wave activity is well organized by outward edges of dense suprathermal ring 
current spots, and the waves are not observed if the ring current ions are not included in 
the real part of dispersion relation. Third, the most intense wave-induced ring current 
precipitation is located in the night MLT sector and caused by modification of the wave 
dispersion relation. The strongest precipitating fluxes of about 8 • 10 6 (cm 2 • s • sr) -1 
are found near L=5.75, MLT=2 during the early recovery phase on 4 May. Finally, the 
nightside precipitation is more intense than the dayside fluxes, even if there are less 



3 


34 intense waves, because the convection field moves ring current ions into the loss cone on 

35 the nightside, but drives them out of the loss cone on the dayside. So convection and 

36 wave scattering reinforce each other in the nightside, but interfere in the dayside sector. 
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37 1 . Introduction 

38 Electromagnetic ion cyclotron (EMIC) waves are a common feature of the Earth 

39 magnetosphere. These waves were observed in the inner [e. g., LaBelle et al. . 1988; 

40 Erlandson and Ukhorskiy, 2001] and outer [Anderson et al, 1992a, b] magnetosphere, 

4 1 at geostationary orbit [Young et al. , 1981; Mauk, 1982], at high latitudes along 

42 the plasmapause [Erlandson et al., 1990], and at ionospheric altitudes [Iyemori and 

43 Hayashi, 1989; Braysy et al, 1998]. Interaction of the ring current (RC) with EMIC 

44 waves causes scattering of ions into the loss cone and leads to decay of the RC 

45 [Cornwall et al., 1970]. This wave-induced RC precipitation was studied widely both 

46 experimentally and theoretically [e. g., Soraas et al, 1999; Erlandson and Ukhorskiy, 

47 2001; Yahnina et al, 2003; Walt and Voss, 2001, 2004; Jordanova et al, 2001; Khazanov 

48 et al, 2002], which produce RC decay times of about one hour or less during the 

49 main phase of storms [Gonzalez et al, 1989]. Obliquely propagating EMIC waves 
so damp due to Landau resonance with thermal plasmaspheric electrons, and cyclotron 
si resonances with thermal, suprathermal, and hot heavy ions [e. .g., Cornwall et al, 

52 1971; Anderson and Fuselier, 1994; Home and Thome, 1997; Thome and Home, 1994; 

53 1997]. Subsequent transport of the dissipating wave energy into the ionosphere causes 

54 ionosphere temperature enhancements [e. g., Gurgiolo et al, 2005]. Cornwall et al. 

55 [1971] employed the mechanism of resonant energy transfer to electrons to explain 

56 stable auroral red arc emissions during the recovery phase of storms. Measurements 

57 taken aboard the Prognoz satellites revealed a “hot zone” near the plasmapause where 
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58 the temperature of core plasma ions can reach tens of thousands of degrees [ Bezrukikh 

59 and Gringauz, 1976; Gringauz, 1983; 1985]. The earliest results regarding the heating 

60 of the cold ions were obtained by Galeev [1975] who considered the induced scattering 

61 of EMIC waves by plasmaspheric protons as an ion heating mechanism. This nonlinear 

62 wave-particle interaction process was used in a plasmasphere— RC interaction model by 

63 Gorbachev et al. [1992], Later, a detailed analysis of thermal ion heating by EMIC 

64 waves was presented by Anderson and Fuselier [1994] and Fuselier and Anderson 

65 [1996]. Relativistic electrons (> 1 MeV) in the outer radiation belt can also interact 
ee with EMIC waves [ Thome and Kennel, 1971; Lyons and Thome, 1972]. Recently, 

67 data from balloon-borne X-ray instruments provided indirect but strong evidence for 

68 EMIC wave-induced precipitation of outer-zone relativistic electrons [Foat et al, 1998; 

69 Lorentzen et al, 2000]. These observations stimulated theoretical and statistical studies 

70 [ Summers and Thorne, 2003; Albert, 2003; Meredith et al, 2003; Loto’aniu et al, 2006] 

71 which demonstrated that EMIC wave-induced pitch-angle diffusion of MeV electrons 

72 can operate in the strong diffusion limit with a time scale of several hours to a day, 

73 and that this mechanism can compete with relativistic electron depletion caused by the 

74 adiabatic effect of Dst during the initial and main phases of a storm. Therefore, EMIC 

75 waves interact well with both the magnet ospheric electrons and ions, and these waves 

76 are strongly influence the particle dynamics in the eV-MeV energy range. 

77 In a number of magnetospheric regimes, a source of free energy for the excitation 

78 of EMIC waves is the temperature anisotropy (T±_ > Tj|) of the hot Fl + distribution 

79 [ Cornwall , 1964, 1965; Kennel and Petschek, 1966], Our understanding of EMIC 
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wave growth and propagation was dramatically changed after measurements on board 
the GEOS 1 and 2 satellites. They revealed the critical role of the thermal He + for 
generation and propagation of EMIC waves [Young et al., 1981; Roux et al, 1982], The 
observations stimulated theoretical studies in which the influence of thermal He + and 
0 + admixtures on EMIC wave properties was considered [Mauk, 1982; Roux et al., 
1982; Rauch and Roux, 1982; Gomberoff and Neira 1983; Gendrin et al, 1984; Denton 
et al. , 1992; Horne and Thorne, 1993]. The effects of energetic RC heavy ions ( He + 
and 0 + ) on the generation of EMIC waves in a multi-ion core plasma (H + , He + , 0 + ) 
were studied by Kozyra et al. [1984]. Horne and Thome [1993] used the “HOTRAY” 
ray tracing program to study the role of propagation and refraction in the generation of 
different branches of EMIC waves in a multi-ion thermal plasma. They found that the 
local growth rate alone cannot determine the resulting wave amplification; propagation 
effects have a major impact on the path-integrated wave gain, and consequently the 
prevalent He + -inode grows preferably at the plasmapause. Recently, Loto’aniu et al. 
[2005] used magnetic and electric field data from the Combined Release and Radiation 
Effects Satellite to obtain the Poynting vector for Pc 1 EMIC waves. They found 
bidirectional wave energy propagation, both away and toward the equator, for events 
observed below 11° |MLAT|, but unidirectional energy propagation away from the 
equator for events outside ±11° of the equator. Engebretson et al. [2005] found a similar 
EMIC wave energy propagation dependence, with mixed direction within approximately 
±20° MLAT, but consistently toward the ionosphere for higher magnetic latitudes. 
These observations allowed Engebretson et al. [2007] to state that “the mixed directions 
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10 2 observed in the above studies near the equator is evidence of wave reflection at the 

103 off-equatorial magnetic latitude corresponding to the ion-ion hybrid frequency. Waves 

104 that reflect would then set up a standing (bi-directional) pattern in the equatorial 

105 magnetosphere. Waves that tunnel through would tend to be absorbed in the ionosphere 
loe and not be able to return to equatorial latitudes.” 

107 Starting from the pioneering work of Kennel and Petschek [1966], it is well-known 
log that the plasma density is one of the most important plasma characteristics controlling 
loo EMIC wave generation; the minimum energy of resonant ions is proportional to the 
no magnetic field energy per particle. In an electron-proton plasma, Cornwall et al. [1970] 
in found that the EMIC wave growth rate maximizes just inside the plasmapause where 

11 2 the Alfven speed is low, falling to zero with both decreasing (because of electron-ion 

113 collisions) and increasing L-shell (because of high critical anisotropy). In the case 
n 4 of a multi-ion magnetosphere, Horne and Thorne [1993] reported a result opposite 
us to that found by Cornwall et al. [1970], namely, the growth rates are substantially 
lie greater outside the plasmapause than just inside the plasmapause. The latter is an 
n 7 effect of heavy ions, and both the above results were reconciled by Kozyra et al. 

ns [1984], However, Horne and Thorne [1993] illustrated that when propagation effects 
ns are properly included, the path-integrated wave gain is indeed larger just inside the 

120 plasmapause. The effect of the plasmapause in EMIC wave generation is very clearly 

121 observed both in experiments [e. g., Fraser and Nguyen, 2001], and in the results of 

122 numerical simulation [Kozyra et al, 1997; Khazanov et al, 2006]. (Of course, the real 

123 magnetospheric situation is more complex, and wave occurrence actually increases with 



124 L shell, which depending on MLT, exhibits a radial structure with a gap between high 

125 and low L-shell events [Anderson et al, 1992 a].) 

126 Recently, Engebretson et al. [ 2007 ] presented the Cluster observations of EMIC 

127 waves in the Pc 1-2 frequency range and associated ion distributions during the October 

128 and November 2003 storms. The most intense waves were observed on 22 November 

1 29 near the end of the rapid recovery phase in the dawn MLT sector at L= 4 . 4 - 4 . 6 . 

130 Generation of these waves was associated with anisotropic RC H + of energies greater 

131 than 10 keV. Although the temperature anisotropy of these energetic protons was high 

132 during the entire 22 November event, EMIC waves were observed only in conjunction 

133 with intensification of the ion fluxes below 1 keV by over an order of magnitude. This 

134 suggests that a suprathermal plasma plays an important role in the destabilization of 

135 the more energetic RC and/or plasma sheet ions, because high energy anisotropic RC 

1 36 and/or plasma sheet proton distributions appeared to be a necessary but not sufficient 

137 condition for the occurrence of EMIC waves. Similarly, studying Pc 1-2 events in the 

1 38 dayside outer magnetosphere, Engebretson et al. [ 2002 ] and Arnoldy et al. [ 2005 ] found 

139 that greatly increased fluxes of low energy protons are crucial for the destabilization of 
mo the anisotropic RC protons. Those observations provide clear evidence that both the 

141 cold plasmaspheric plasma (and, of course, heavy ion content) and the suprathermal 

142 1 keV) ions injected from the plasma sheet (and/or ion outflow from the ionosphere) 

143 control EMIC wave excitation in the RC. On the other hand, an assumption that the 

144 total plasma density /composition is dominated by the thermal plasma was made in 

145 previous RC-EMIC wave modeling efforts, and RC ions were not included in the real 
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we part of the wave dispersion relation [Kozyra et al., 1997; Jordanova et al., 1998b, 2001; 
H7 Khazanov et al, 2006], but only in the EMIC wave growth rate. As a result, EMIC waves 
ms are only generated near the plasmapause in all these theoretical models. Consequently 

149 we generalize our previous self-consistent RC-EMIC wave model [Khazanov et al., 2006] 

150 to take into account the effect of RC ions in the real part of the EMIC wave dispersion 

151 relation. 

152 The present study further develops a self-consistent theoretical model of RC and 

153 propagating EMIC waves in a multi-ion magnetospheric plasma [Khazanov et al., 2006], 

154 where we take into account the RC ions in the real part of dispersion relation for the 

155 He + - mode. This article is organized as follows: In section 2 we provide the system of 
us equations which govern our global theoretical model, as well as the initial/boundary 

157 conditions used in the simulation of the May 1998 storm; In section 3 we present both 

158 the spatial distribution of the total plasma density (thermal + higher energies) during 

159 the May 1998 event, and the fine energy structure of the RC phase space distribution 
iso functions; In section 4, the effect of plasma density on the EMIC wave growth is 

161 illustrated; In section 5, role of the RC ion thermal effects in the He + - mode dispersion 

162 relation is analyzed; In section 6, results of simulation are presented; Finally, in section 7 

163 we summarize the new features of the model, and the findings of the paper. 
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164 2. Equations of Global Model, Approaches and 
i66 Initial/Boundary Conditions 


For RC species H + , 0 + . and He + , we simulate the RC dynamics by solving the 
bounce-averaged kinetic equation for the phase space distribution function (PSDF), 
F(ro, <p, E, no, t). The PSDF depends on the radial distance in the magnetic equatorial 
plane r 0 , geomagnetic east longitude <p, kinetic energy E, cosine of the equatorial pitch 
angle p 0 , and time t [see, e. g., Fok et al, 1993; Jordanova et al., 1996]. We use the 
bounce-averaged kinetic equation for the H e + -mode of EMIC waves to describe the 
wave power spectral density. This equation was originally derived by Khazanov et al. 
[2006], and explicitly includes the EMIC wave propagation, refraction and reflection in 
a multi-ion magnetospheric plasma. Following to Khazanov et al. [2006], we ignore the 
slow azimuthal and radial drifts of the waves during propagation, and use the reduced 
wave kinetic equation. So the resulting system of governing equations take the form: 
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lee In the left-hand side of equation (1), all the bounce-averaged drift velocities are denoted 
i 67 as {•••), and may be found in previous studies [Jordanova et al., 1994; Khazanov et al. , 
las 2003]. In equation (2), u and 9 0 are the wave frequency and equatorial wave normal 
i 69 angle, respectively, { 9 0 ) is the bounce-averaged drift velocity of the equatorial wave 



11 


170 normal angle, B w is the EMIC wave magnetic field, and (7} is a result of averaging 

171 of the local growth/damping rates, which includes both the wave energy source due 

172 to interaction with RC ions and the energy sink due to absorption by thermal and 

173 hot plasmas, along the ray phase trajectory over the wave bounce period. Note that 

174 equation (2) is accompanied by a system of the ray tracing equations which are not 

175 written here (for details see Khazanov et al. [ 2006 ] and references therein). 

176 The term in the right-hand side of equation ( 1 ) includes losses from charge 

177 exchange, Coulomb collisions, ion-wave scattering, and precipitation at low altitudes 
ns [Jordanova et al., 1996 , 1997 ; Khazanov et al., 2002 , 2003 ]. Loss through the dayside 
179 magnetopause is taken into account allowing a free outflow of the RC ions from a 

wo simulation domain. The bounce-averaged pitch angle diffusion term in the right-hand 

181 side of equation ( 1 ) is a functional of the EMIC wave power spectral density, B^, i. e. 

182 the diffusion coefficient has the form = (D M0;Al0 (B^(-))). On the other hand, 

1 83 (7) in equation ( 2 ) is a functional of the phase space distribution function, F, i. e. 

1 84 (7) = (7 (F (•))}. So equations ( 1 ) and ( 2 ) self-consistent ly describes the interacting 

185 RC and EMIC waves in a quasilinear approximation. It should be emphasized that in 
iso order to describe the wave-particle interaction in equation (1) we have to know the 

i 87 off-equatorial power spectral density distribution for EMIC waves, and this distribution 
las can then be mapped from the magnetic equator using solutions of the ray tracing 

189 equations. 

190 The geomagnetic field in our simulation is taken to be a dipole field. The electric 

191 field is expressed as the shielded (exponent 2 ) Volland-Stern convection field [ Volland, 
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W 2 1973; Stern, 1975] which is Ap-dependent, with a corotation field [see, e. g., Lyons and 

193 Williams, 1984]. The equatorial thermal electron density distribution is calculated with 

194 the time-dependent model of Rasmussen et al. [1993]. For modeling the RC-EMIC 

195 wave interaction and wave propagation we also need to know the density distribution 
we in the meridional plane, in the present study we employ an analytical density model 
197 which includes the product of three terms; ( 1 ) diffusive equilibrium model term 

os [ Angerami and Thomas, 1964], (2) lower ionosphere term, and (3) plasmapause and 

199 outer magnetosphere term. This analytical model is adjusted to the Rasmussen model 

200 at the equator. So the resulting plasmaspheric density model provides a 3D spatial 

201 distribution for electrons, and an ion content assumed to be 77% for H + , 20% for He + , 

202 and 3% for 0 + . Geocoronal neutral hydrogen number density, needed to calculate 

203 loss due to charge exchange, is obtained from the spherically symmetric model of 

204 Chamberlain [1963] with its parameters given by Rairden et al. [1986]. 

205 In order to study Dst variation during the May 1998 storm period, and to calculate 

206 the energy content for the major RC ion species, H + , 0 + , He + , Farrugia et al. [2003] 

207 used the RC kinetic model of Jordanova et al. [1998a]. They found that during this 

208 storm the energy density of H + is greater than twice that of 0 + at all MLTs, and 

209 the contribution of He + to the RC energy content is negligible. This implies that 

210 RC 0 + content do not exceed 30% during the main phase of this storm. Note that 
2 u above estimation was obtained from a simulation without oxygen band waves. On the 

212 other hand, Brdysy et al. [1998] observed very asymmetric 0 + RC during the main 

213 phase of the April 2-8, 1993 storm, which suggests that the RC oxygen ion loss rate is 
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considerably faster than the drift speed. This result is difficult to explain in terms of 
charge exchange and Coulomb scattering, and suggests that the production of EMIC 
waves contributes significantly to RC 0 + decay during the main and early recovery 
phases. In other words, due to generation of the 0 + -mode EMIC waves, most RC 0 + 
precipitates before reaching the dusk MLT sector [Braysy et al., 1998], Therefore, to 
estimate the RC 0 + content correctly, the 0 + -mode should be included in simulation, 
and it is likely that Farrugia et al. [2003] overestimated the RC 0 + content during 
May 1998. Anyhow, the calculations of Thorne and Home [1997] clearly confirm that 
the above RC 0 + percentage cannot significantly suppress ATe + -mode amplification, 
and only slightly influences the resulting wave growth. It is for this reason we chose to 
initially exclude RC 0 + in our particular simulation of May 2-7, 1998, and to assume 
that the RC is entirely made up of energetic protons. 

The night-side boundary condition is imposed at the geostationary distance in 
our model, and we use the flux measurements during the modeled event obtained from 
the Magnetospheric Plasma Analyzer and the Synchronous Orbit Particle Analyzer 
instruments on the geosynchronous LANL satellites. Then, according to Young et al. 
[1982], we divide the total flux measured at geostationary orbit between the RC II + , 
0 + , and He + depending on geomagnetic and solar activity as measured by Kp and F l07 
indices. Only the H + fluxes were used as a boundary condition in the simulation. 

To obtain the self-consistent initial conditions for equations (1) and (2), the 
simulation was started at 0000 UT on 1 May, 1998 using a background noise level for 
the He + - mode of EMIC waves [e. g., Akhiezer et al., 1975b], the statistically derived 
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236 quiet time RC proton energy distribution of Sheldon and Hamilton [1993], and the 

237 initial pitch angle characteristics of Garcia and Spjeldvik [1985]. The initial the RC 

238 and EMIC wave distributions are derived independently, and of course, have nothing 

239 to do with a particular state of the magnetosphere during a simulated event. Only 

240 the boundary conditions provided by the LANL satellites can be considered as data 

241 reflecting a particular geomagnetic situation (and, to a certain extent, the employed 

242 plasmasphere and electric field models driven by Kp). Therefore, before simulation of 

243 a particular geomagnetic event can be possible, we first seek an initial state for the 

244 RC and EMIC waves that is self-consistent and reflects the particular geomagnetic 

245 situation. In our case, this was done by running the model code for 24 hours. In about 

246 20 hours of evolution, the wave magnetic energy distribution reaches a quasistationary 
24 ? state indicating that the RC-EMIC wave system achieves a quasi-self-consistent state. 

248 (Note that 20 hours has nothing to do with the typical time for wave amplification and 

249 instead reflects the minimum time needed to adjust RC and waves to each other and 

250 to the real prehistory of a storm.) So the self-consistent modeling of the May 1998 

25 1 storm period is started at 0000 UT on 2 May (24 hours after 1 May 0000 UT) using 

252 solutions of equations (1) and (2) at 2400 UT on 1 May as the initial conditions for 

253 further simulation. 



15 


254 


255 


256 


257 


258 


259 


260 


261 


262 


263 


264 


265 


266 


267 


268 


269 


270 


271 


272 


273 


274 


3. Distribution of Plasma Density and Energy Structure of RC 
PSDFs 


3.1. Spatial Patterns of Plasma Density During the May 1998 Storm 


From the results of our simulation we select seven snapshots which represent the 
intervals of the most enhanced plasma sheet H + injection into the RC region. The 
selected equatorial plasma density distributions are presented in Figure 1. The first row 
in this Figure shows the electron plasma density distribution from the Rasmussen et al. 
[1993] model, and the second row provides a sum of the corresponding plasma density 
from the first row and the RC H + density. Note that starting from high L-shell, the RC 
ions dominate the thermal plasma excepting a plasmaspheric drainage plume, and below 
we shell concentrate only on cases of pronounced density enhancement during plasma 
sheet ion injections. The first plasma sheet ion injection appears about 32 hours after 1 
May, 0000 UT (not shown), which affects the density distribution for about 16 hours, 
while the RC ions only slightly modify the plasma density distribution after 48 hours 
(not shown). During this interval, the RC H + density dominates the thermal plasma 
in the dusk-midnight MLT sector (see hours 33 and 34 in Figure 1). The second ion 
injection starts about 56 hours (not shown). The snapshots at hour 60 show the most 
distinct pattern of the cold and total plasma density during this injection event when 
the RC H + dominates the thermal plasma density in the nightside through the entire 
dusk-dawn MLT sector. Again, there are only minor differences between the density 
snapshots at 68 hours (not shown). The third plasma sheet ion injection shown in 


Figure 1 
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275 Figure 1 starts at about 76 hours and impacts the plasma density distribution through 

276 hour 90 (not shown). This injection is most intense comparing to previous ones, and 

277 the RC H + dominance is observed in the greatest L-shell and MLT extents encircling 

278 a great part of the globe during the third injection. The results of our simulation are 

279 in qualitative agreement with the RC density distribution obtained by Zaharia et al. 

280 [2006] during the moderate geomagnetic storm of 21-23 April 2001. 

28 1 We presented only the RC H + density distribution above, and did not say 

282 anything about the distribution of the electron density. It is obvious that in all “slow” 

283 magnetospheric processes the quasi-neutrality condition should hold. This implies 

284 that electrons have the same density distribution as the ions. Quasi-neutrality can 

285 be sustained by both the energetic plasma sheet electrons injected along with ions, 

286 and/or the cold ionospheric electrons due to field-aligned currents. The resulting 

287 electron temperature strongly affects the Coulomb energy degradation of the RC ions, 

288 the resonant Landau damping of EMIC waves, and barely influences the EMIC wave 

289 dispersive properties (see, e. g., Khazanov et al. [2007], Akhiezer et al. [1975a]). 

290 Khazanov et al. [2007] demonstrated that both the EMIC wave Landau damping and 

291 collisional RC energy dissipation are maximized for an electron temperature about 

292 1 eV. This is the temperature adopted in our RC-EMIC wave model for thermal plasma 

293 [Khazanov et al. , 2003], Therefore, if we do not track the electron dynamics and keep 

294 T e — 1 eV for the entire simulation domain, we can potentially underestimate the EMIC 

295 wave energy, especially at high L-shells during the main and recovery storm phases when 

296 RC ions dominate the thermal plasma. Below we assume that plasma is quasi-neutral 
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and that the electron temperature is 1 eV throughout the entire simulation domain 
during the May 1998 event. 


3.2. Fine Energy Structure of RC PSDFs 


The new RC ions, injected from the plasma sheet in the night MLT sector, cause 
impressive plasma density enhancement for high L-shells during the main and recovery 
storm phases. This feature is clearly observed in our simulation, but in Figure 1 we 
presented only the RC H + density distribution, and did not analyze the fine PSDF 
energy structure. To consider the energy distributions of the RC H + , we selected four 
representative cases among the snapshots in Figure 1. The corresponding PSDFs are 
shown in Figure 2. All the PSDFs are taken in the equatorial plane, and integrated over 
the entire solid angle, while the effective RC proton temperature parallel to geomagnetic 
field line, T\\, is calculated for the entire energy range (100 eV — 430 keV). In order to 
more clearly demonstrate change in the PSDF slope, we use a linear energy scale in a 
low energy domain of the distribution, whereas the high energy part is depicted with 
a logarithmic energy scale. As follows from the left-hand side of Figure 2, there is a 
transition region in all the PSDFs which separates relatively warm ions from the more 
hot and tenuous component. (The transition from a steep profile to more horizontal 
profile corresponds to the transition from a small to a higher effective ion temperature.) 
So we observe at least two ion populations which constitute the plotted RC ion PSDFs; 
(1) the dense and relatively cold low energy RC component, and (2) the rare and 
hotter high energy RC component. The boundary between these two ion components 


Figure 2 
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318 is located at slightly different energy depending on each case, which from Figure 2, is 

319 about 1 — 1.5 keV. Note that PSDFs at hours 80 and 82 include, respectively, four and 

320 three ion populations with different effective temperatures; the PSDF taken at hour 80 

321 changes slope at energies near 1, 10, and 130 keV, whereas the PSDF at hour 82 changes 

322 slope near 0.5 and 20 keV. So the results in Figure 2 clearly demonstrate that plasma 

323 density modification due to the plasma sheet H + injection into the RC region is mostly 

324 caused by low energy ions with energy ^ 1 keV. 


325 4 . Effect of Plasma Density on EMIC Wave Growth 

326 The effective proton temperatures transverse to Tj_, and along Xy, the geomagnetic 

327 field line, comply with the inequality T± > Tj| in many space plasma regimes. If the ion 

328 temperature anisotropy, A = T±/T\\ — 1, exceeds some positive threshold, EMIC waves 

329 can be unstable [ Kennel and Petschek, 1966; Cornwall et al, 1970]. The growth rate 

330 for these waves critically depends on the characteristic energy for cyclotron interaction, 

331 which, as defined by Kennel and Petschek [1966], is just the local geomagnetic field 

332 energy per particle, having the form E c = B 2 / ( 8 ivn e ). So, according to Kennel and 

333 Petschek [1966], the local growth rate for EMIC waves should be particularly sensitive 

334 to the local plasma density. Assuming that the RC is entirely made up of energetic H + , 

335 Figure 3 plots the dependence on plasma density of the local equatorial growth/damping 

336 rate for the He + mode EMIC waves. Note that the calculated growth/ damping rates 

337 in Figure 3 are due to the RC wave interaction only, and the wave absorption due to 

338 thermal plasma is omitted (but, of course, this effect is included in global simulation). 


Figure 3 
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339 All the results in Figure 3 are obtained for the wave frequency v == 0.475 Hz, and case 

340 (a) is just taken from our global model without any modification at location L=5.25, 

341 MLT=15 at 48 hours (n e — no = 68.3 cm -3 , and B = 215.3 nT). In order to produce 

342 the results (b), (c), and (d), we need only re-normalize the local plasma density as 

343 n e = 1.2 x n 0 , n e — 1.5 x n 0 , and n e = 2.0 x no, respectively. As follows from Figure 3, 

344 transitioning from case (a) to case (b) increases the peak growth rate by a factor 1.4, 

345 extends the region of growth, and makes the wave damping negligible. Further increase 

346 of the number density eliminates the region of wave damping. According to [Kennel and 

347 Petschek, 1966], the growth rate dependence on plasma density is 7 ~ exp (— l/n e )/ v /?v 

348 So, although the characteristic energy decreases with increasing plasma density, the 

349 growth rate can both increase or decrease depending on the wave normal angle (see 

350 Figure 3). For a particular wave normal angle, it depends on whether we move to the 

351 growth rate maximum with density increase or whether we move from the maximum. 

352 5 . Effects of RC Temperature on EMIC Wave He + — Mode 

Although the results presented in subsection 3.2 clearly demonstrate that the 
observed plasma density enhancement is caused by a low energy (,<; 1 keV) population 
of the RC, this does not allow us to evaluate the effects of the RC ion temperature on 
the EMIC wave dispersive properties. In order to characterize the temperature effects 
in the EMIC wave dispersion relation, we use the following parameters [see, e. g., Stix. 
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1992; Akhiezer et al., 1975a] 


A* 


( k±v±,i 

V V2 a 


, Ci 


to dh fi* 
k\\v u 


, i = e : H + ,He + ,0 + , 


( 3 ) 


where fij is the particle gyrofrequency, and k±_ (uj.,* = y 2Xj_ i j/m i ) and fty 


27]| (i /mi) are the components of the wave normal vector (thermal velocity) transverse 
to and along geomagnetic field lines, respectively; A * is the squared ratio of Larmor 
radius to transverse wave length; and Q is the squared ratio of longitudinal wave length 
to a typical particle displacement along the field line during a wave period. The finite 
Larmor radius effects are negligible if A * << 1. On the other hand, the plasma particles 
become unmagnetized if A, >> 1, and as a consequence the external magnetic field 
disappears in the wave dispersion relation. So the Larmor radius effects are most 
important for an intermediate case when the wave and particle parameters give A* ~ 1. 
The magnitude of Q not only characterizes the importance of “longitudinal” thermal 
effects, but also determines the effectiveness of the resonant wave damping/ growth. For 
instance, the number of resonating particles is small if Q » 1, and as a result, plasma 
waves can exist for a long time without substantial damping. So the role of thermal 
effects in the wave dispersion relation depends on the magnitude of both Q and A,. For 
example, if these parameters comply with the inequalities A, << 1 and Q >> 1, in many 
cases (but not always!) the leading term in a real part of dispersion relation still comes 
from a cold plasma approximation (limit A* = 0 and Q — » oo, e. g., Stix [1992]). So 
depending on the magnitudes of Q and A *, the thermal terms may be a minor correction 
only, or they can dominate the “cold plasma limit” term. 
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Until now, we discussed only the RC H + . Although the RC H 1 dominate both 0 + 
and He + during the May 1998 storm [Farrugia et al. 2003], and we do not simulate the 
RC 0 + and He + in the present study, the heavy ions participate in the RC dynamics 
and can influence the magnetospheric heavy ion content, especially during the main and 
early recovery storm phases. Despite the importance of the hot heavy ions for the EMIC 
wave characteristics (see, e. g., Kozyra et al [1984]), in all previous studies we assumed 
that the total ion composition is dominated by the ion composition of the thermal 
plasma and did not take into account the RC ions in the real part of the wave dispersion 
relation [Khazanov et al, 2002, 2003, 2006, 2007], including the RC ions in the imaginary 
part only. In all those papers, when we described the EMIC wave dispersive properties 
we used the electron density distribution from the time-dependent Rasmussen et al. 
[1993] model, and the ion content was assumed to be 77% for H + , 20% for He + , and 
3% for 0 + . (Although the assumed ion content is in the range of 10 — 30% for He + 
and 1 — 5% for 0 + following observations by Young et al [1983, 1977] and Horwitz et 
al [1981], it only approximately describes the real ion percentage and, of course, does 
not reflect its variability, especially during the magnetically active periods.) Now we are 
going to take into account the RC ions in the real part of the EMIC wave dispersion 
relation which can strongly modify the heavy ion percentage. In spite of this, for the 
purpose of comparison with previous results, we keep the earlier adopted ion percentage 
(77% for H + , 20% for He + , and 3% for 0 + ) throughout the entire simulation domain 
even if this percentage is mainly determined by the suprathermal/hot ion composition. 

It follows from equation (3), assuming that all the RC ions ( H + , He + , 0 + ) have 
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nearly the. same temperature, that parameters A* relate to each other as masses of 
the corresponding RC ions. Then, considering the most dense suprathermal spots in 
Figure 1, we find that for the He + - mode the following inequality 

Xff+ < Xff e + < A o+ « 1 << C,He+ rS Co+ << C H+ (4) 

holds. Note that in order to obtain inequalities (4), we used vx,% and v\\ ti calculated for 
the entire energy range; parameters A $ and Q could be even closer to the cold plasma 
limit if all the effective temperatures are calculated for a low energy RC component 
only (see subsection 3.2), which gives the greatest contribution to the plasma density 
enhancement observed in night side during the main and recovery storm phases. In the 
limit (4), the structure of thermal terms in the EMIC wave dispersion equation can be 
found, e. g., in [Stix, 1992; Akhiezer et al, 1975a] where the finite Larmor radius effects 
may be omitted. The greatest thermal term (A||) in the dispersion equation for the 
EMIC wave H e + -mode comes from the RC H + during the May 1998 storm with the 
following ranking 

A jl Cff + )»A,(0+)~A||(.ffe+). (5) 

393 So only term A\\(H + ) can potentially compete with the “cold plasma limit” term in 

394 the He + - mode dispersion equation. Considering the most dense suprathermal spots 

395 in Figure 1 , we find that A||(Ef + ), as a rule, can be neglected in comparison with the 

396 “cold” term in the i?e + -mode dispersion relation. 
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6. Results and Discussions 

Summarizing all the assumptions and conclusions we did in sections 3 and 5: 

(1) Plasma is quasi-neutral (see subsection 3.1); (2) the electron temperature is 
1 eV through the entire simulation domain (subsection 3.1); (3) the plasma density 
enhancement observed in Figure 1 is caused by a low energy (^ 1 keV) population of 
the RC ions (subsection 3.2), while the RC H + ions dominate both the RC 0 + and 
He + during May 1998; (4) the ion percentage is 77% for H + , 20% for He + , and 3% 
for 0 + through the entire simulation domain (section 5); and (5) the thermal effects of 
electrons and the RC ions may be neglected in the real part of the 77 e + -mode dispersion 
relation (see subsections 3.1 and section 5). 

6.1. Global Distribution of He + — Mode 

The equatorial (MLT, L-shell) distributions of the squared wave magnetic field, 

fWmax / >7r _ 

Bl(r 0 ,(p,t)= du dOoB^ (r 0 , <p, t, u, 6 q) , (6) 

^min ''O 

are shown in Figure 4 for the He + - mode of EMIC waves. These simulation results 
are based on the system of governing equations (1) and (2) along with the ray tracing 
equations. The results in the first row are obtained when the RC ions are only treated 
as a source of free energy to generate EMIC waves, and omitted in the real part of the 
wave dispersion relation. The second row shows the case when the RC ions are taken 
into account in both the real and imaginary parts of the wave dispersion relation. There 
is an essential difference between the EMIC wave energy distributions in the first and 


Figure 4 
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second rows. Modification of the EMIC wave dispersive properties due to RC ions leads 
to a relatively minor spatial redistribution of the “old” wave active zones presented in 
the first row, and mainly alters the wave intensities. The qualitative difference between 
the first and second rows appears during the recovery phase in the postmidnight-dawn 
MLT sector for L > 4.75 (hours 82 and 84). In these regions, “new” EMIC waves are 
generated due to modification of the wave dispersion by RC, and we do not observe any 
wave activity in corresponding snapshots in the first row. The B-field distributions are 
organized by the locations of sharp gradient in the total density of thermal plasma and 
RC as expected from previous studies [Horne and Thorne, 1993; Khazanov et al, 2006]. 
(The sharp density drop counteracts the refraction caused by the magnetic field gradient 
and curvature. As a result, net refraction is suppressed, and the ife + -mode grows 
preferentially at these locations.) At the same time, we note that a radial extension of 
wave zones in the second row is slightly greater than that in the first row. 

Let us now discuss the new feature caused by the modified EMIC wave dispersion 
and clearly observed in Figure 4. Recently, Engebretson et al. [2007] presented 
measurements of EMIC waves in the Pc 1-2 frequency range and the associated ion 
distributions obtained Cluster. During the October and November 2003 magnetic 
storms, the most intense waves were observed on 22 November near the end of a rapid 
recovery phase from 0825 to 0850 UT; located near dawn for L=4.4-4.6 and at an 
average ML AT 18°. The waves were primarily transverse, propagated away from 
the equator, and predominantly left-hand polarized. Compared to the local proton 
gyrofrequency, these waves had a normalized frequency of X=0.34, somewhat higher 
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than the local He + gyrofrequency (X=0.25). The free energy to generate those waves 
was associated with anisotropic RC IT r of energies greater than 10 keV. Note that the 
upper energy range of increased energy fluxes may well extend beyond the 40 keV limit 
of the Cluster CIS instrument. Although the temperature anisotropy of these energetic 
(> 10 keV) protons was high during the entire 22 November pass, EMIC waves were 
observed only in conjunction with intensification of the ion fluxes below 1 keV by over 
an order of magnitude. This suggests that the suprathermal plasma plays an important 
role in the destabilization of the more energetic RC and/or plasma sheet ions, and the 
high energy anisotropic RC and/or plasma sheet proton distributions appeared to be 
a necessary but not sufficient condition for the occurrence of EMIC waves. Similarly, 
studying Pc 1-2 events on the dayside outer magnetosphere, Engebretson et al. [2002] 
and Arnoldy et al. [2005] found that greatly increased fluxes of low energy protons are 
crucial for the destabilization of the high energy anisotropic RC protons. 

The satellite observations by Engebretson et al. [2007] support our theoretical 
results presented in Figure 4. Indeed, in the second row we see intense EMIC waves (up 
to a few nT 2 ) in the postmidnight-dawn sector (for L > 4.75) during the recovery phase 
from 82 to 84 hours. This wave activity is not observed if the RC ions are not included 
in the real part of the wave dispersion relation (compare the first and second rows in 
Figure 4). At the same time, we note that Engebretson et al. [2007] observed waves with 
a normalized frequency X=0.34, whereas we consider the He + - mode of EMIC waves 
with X < 0.25. (The most intense burst of Pc 1 waves studied by Arnoldy et al. [2005] 
was measured by the Polar satellite with a local normalized frequency of X=0.2, so the 
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459 waves were also He + - mode.) For the purpose of comparison with previous results, in 
46d the present study we kept the ion percentage the same as in our earlier studies, namely, 

461 77% for H + , 20% for He + , and 3% for 0 + . Then the most effective generation takes 

462 place for the He + - mode in the frequency range Qo+ < u) < 0/i e + [see, e. g., Kozyra et 

463 al, 1984; Horne and Thome, 1993; Khazanov et al, 2003], (Note that only waves in 

464 the left-hand polarized part of the dispersive surface can grow, and the corresponding 

465 wave frequencies should be in the range between the cross-over frequency and Qfj e +.) 

466 This heavy ion content, however, differs strongly from the ion percentage reported by 
46? Engebretson et al. [2007], For example, they observed 81% of H + , 3% of II e + , and 16% 

468 of 0 + on November 22, 2003 at 0740 UT, qualitatively different from the percentage 

469 we used in the simulation. Such a great amount of RC 0 + , in combination with small 

470 amounts of He + , should suppress the He + - mode, and conversely favor the H + -mode. 

471 Self-consistent modeling of the /7 + -mode is beyond the scope of the current study, and 

472 should be done separately. (Strictly speaking, EMIC waves are very sensitive to the 

473 the heavy ions, so wave simulation requires more realistic dynamic models of the global 

474 distribution for each ion species which, unfortunately, are currently not available.) At 

475 present, we believe that the crucial role of low energy RC and/or plasma sheet protons 

476 in the destabilization of the high energy anisotropic RC protons is well established both 

477 experimentally and theoretically. We also think that this feature depends on the wave 

478 mode only quantitatively, and the qualitative effect itself does not depend on the wave 

479 mode. 
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480 6.2. Wave-Induced RC Precipitation 

One of the most pronounced consequences of the RC-EMIC wave interaction is 
the scattering of RC ions into the loss cone. This process is one of the processes that 
lead to decay of RC [see, e. g., Cornwall et al., 1970], especially during the main and 
early recovery phases of storms when decay time of about one hour or less is possible 
[Gonzalez et al. , 1989]. The EMIC wave-induced RC precipitation was studied widely 
both experimentally and theoretically [e. g., Erlandson and Ukhorskiy, 2001; Yahnina 
et al. , 2003; Walt and Voss, 2001, 2004; Jordanova et al, 2001]. Although the effect 
of EMIC waves on RC ion precipitation during the May 1998 storm was discussed 
previously [e. g., Khazanov et al., 2002, 2007], we present a few precipitating patterns 
that demonstrate the new features caused by modification of the EMIC wave dispersion 
relation. The RC precipitating flux is calculated as 

1 rE 2 r\ /*1 

Jlc= rT lc = / d Ah), (7) 

481 where q/ c is the cosine of the equatorial pitch angle at the boundary of loss cone, and 

482 j is the equatorial ion differential flux. In Figure 5 we show selected snapshots of the 

483 precipitating fluxes integrated over the energy range 1 — 50 keV. As before, the first 

484 row shows the results without the RC ions in the real part of the EMIC wave dispersion 

485 relation, while the second row shows precipitation when the RC ions are taken into 

486 account in both the real and imaginary parts of the wave dispersion relation. There are 

487 many differences between the first and second rows. The most intense ion precipitation 

488 is due to “new” wave activity, and located in the night MLT sector. The strongest 


Figure 5 




28 


489 


490 


491 


492 


493 


494 


495 


496 


497 


498 


499 


500 


501 


502 


503 


504 


505 


506 


507 


508 


509 


fluxes of about 8 • 10 6 (cm 2 • s • sr) -1 are observed near L=5.75, MLT=2 during the early 
recovery phase of the storm (see hour 82 in Figure 5). This precipitation is two times 
greater than a greatest flux from a previous study of the May 1998 storm by Khazanov 
et al [2007]. The very interesting result can be derived by comparing Figure 5 with 
Figure 4; the wave-induced night side precipitation is more intense than the day side 
fluxes, even if there are less intense waves (compare locations L=4.5, MLT=16, and 
L=5.75, MLT=2 in the 82 hour snapshots). The major reason for this feature is a 
magnetospheric convection field which acts oppositely in day and in night sides moving 
RC ions into the loss cone on the nightside, and driving them out of the loss cone 
on the dayside. So the magnetospheric convection and the wave scattering reinforce 
each other on the nightside, but subtract on the dayside. Of course, we have to recall 
that characteristics of the wave normal angle distribution can strongly impact the 
effectiveness of RC ion scattering [ Khazanov et al, 2007]. 

7. Conclusions 

In this paper we have further developed a self-consistent model of RC ions and 
propagating EMIC waves by Khazanov et al. [2006], We have taken into account RC 
ions in the real part of dispersion relation for the i7e + -mode of EMIC waves. This is a 
new feature of the present model and generalizes the limiting assumption that the total 
plasma density was dominated by the thermal plasma made by all previous RC-EMIC 
wave models, so that the RC ions were not taken into account in the real part of the 
wave dispersion relation [ Kozyra et al, 1997; Jordanova et al, 1998b, 2001; Khazanov 
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5io et al, 2003, 2006] but only in the imaginary part, i. e., in the EMIC wave growth rate. 
5u This assumption is not always valid, especially for high L-shells during the main and 

512 recovery storm phase when the newly injected RC ions dominate the thermal plasma 

513 (see results of our simulation in Figure 1). Recent satellite observations during the 

514 November 2003 magnetic storm by Engebretson et al. [2007] showed that although 

515 the temperature anisotropy of energetic (> 10 keV) RC protons was high during the 

516 entire 22 November 2003 perigee pass, EMIC waves were observed only in conjunction 

517 with intensification of the ion fluxes below 1 keV by over an order of magnitude. This 
sis suggests that the suprathermal plasma (^<C 1 keV) plays an important role in the 

519 destabilization of the more energetic RC and/or plasma sheet ions such that high energy 

520 anisotropic RC and/or plasma sheet proton distributions appeared to be a necessary 

521 but not sufficient condition for occurrence of EMIC waves. 

522 To demonstrate the role of RC ions in the real part of EMIC wave dispersion 

523 relation, we have simulated the May 1998 storm, and have presented and discussed 

524 the global distributions of the total plasma density, the energy of the He + ~ mode, and 

525 the wave-induced RC precipitation. The main conclusions of our simulation can be 

526 summarized as follows. 

527 1. The new RC ions, injected from the plasma sheet in the night MLT sector, causes 

528 plasma density enhancements for high L-shells during the main and recovery storm 

529 phases. This feature is clearly observed in our simulation (see Figure 1), and the plasma. 

530 density enhancement is mostly caused by the suprathermal H + (rS 1 keV). 


531 


2 . During the recovery phase, modification of the wave dispersion relation by RC 
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ions leads to a dramatic change in 'the wave patterns in the nightside MLT sector for 
L > 4.75. 

3. The Cluster observations of EMIC waves and associated ion distributions during 
the November 2003 magnetic storm [Engebretson et al. , 2007] support our theoretical 
results presented in Figure 4. In the second row of Figure 4 we see intense EMIC waves 
(up to a few nT 2 ) in the postmidnight-dawn sector during the recovery storm phase 
from 82 to 84 hours. This wave activity is not observed if the RC ions axe not included 
in the real part of the wave dispersion relation (compare the first and second rows in 
Figure 4). 

4. The most intense wave-induced RC precipitation is due to modification of the 
wave dispersion relation, located in the night MLT sector. The strongest precipitating 
fluxes of about 8 • 10 6 (cm 2 • s • sr) -1 are observed near L=5.75, MLT=2 during the early 
recovery phase of the storm (see hour 82 in Figure 5). The wave-induced nightside 
precipitation is more intense than the dayside fluxes, even if there are less intense 
waves (compare the results at L=4.5, MLT=16, and L=5.75, MLT =2 in the 82 hour 
snapshots) . 
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Figure 1. Equatorial plasma density distributions during the May 1998 event. The 
first row shows the cold electron plasma density distribution from the Rasmussen et al. 
[1993] model, and the second row provides a sum of cold plasma density and RC H + 

790 

density as it follows from the simulation. The first, the second, and the third plasma 

sheet ion injections affect the total density distribution during 33-48, 58-68, and 78-90 

hours, respectively. The specified hours are counted from 0000 UT on 1 May, 1998. 
Figure 2. Simulated phase space distribution function for the RC H + . All the PSDFs 

are shown in the equatorial plane, and integrated over the entire solid angle. For each 

PSDFs, the first and the second numbers in parenthesis are the L-shell and MLT location, 

791 

respectively. The corresponding RC proton temperature along the geomagnetic field line, 
Tj] , is calculated for the entire energy range. Note that there are the linear and logarithmic 
energy scales in the left-hand and right-hand boxes, respectively. 

Figure 3. Equatorial growth/ damping rates versus the wave normal angle for the He + - 
mode of EMIC waves. The RC is assumed to be entirely made up of energetic protons, 
the thermal plasma consists of the cold electrons, and 77% of H + , 20% of He + , and 3% 
of 0 + . and the wave resonate interaction with thermal plasma is omitted. All the results 
are obtained for the wave frequency u = u>/ 2 tt — 0.475 Hz, and taken from our global 

792 

model at location L=5.25, MLT=15 (£> = 215.3 nT), at 48 hours after 1 May 1998, 
0000 UT. (a) The electron number density is also determined by the global model, and 
n e — no = 68.3 cm -3 (nominal case). In order to produce the results (b), (c), and (d), 
we keep all parameters the same, except the electron number densities n e — 1.2 x n ih 
n e = 1.5 x n 0 , and n e = 2 x n 0 are respectively adopted. 
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Figure 4. Snapshots of the equatorial (MLT, L-shell) distributions of squared wave 
magnetic field for the /TtCmode. The results are obtained by solving equations (1) and 
(2) along with the ray tracing equations. The first row corresponds to the case when the 
RC ions are only treated as a source of free energy to generate waves, and omitted in 

793 

the real part of the wave dispersion relation. The second row demonstrates distribution 

when the RC ions are taken into account in both the real and imaginary parts of the 

wave dispersion relation. In both cases, the total ion composition is assumed to be 77% 

of H + , 20% of He + , and 3% of 0 + through an entire simulation domain. 

Figure 5. The RC proton precipitating fluxes averaged over the equatorial pitch-angle 

loss cone and integrated over the energy range 1 — 50 keV. The first row represents the 

794 results without the RC ions in the real part of the EMIC wave dispersion relation. The 
second row shows precipitation in a case when the RC ions are taken into account in 
both the real and imaginary parts of the wave dispersion relation. 
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